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We consider the Asakura-Oosawa-Vrij model of hard sphere colloids and ideal polymer coils in 
contact with a planar hard wall at (colloidal) liquid-gas coexistence. Using extensive numerical 
density functional calculations, the liquid-gas, wall-liquid and wall-gas interfacial free energies are 
calculated. The results are inserted into Young's equation to obtain the contact angle between the 
liquid-gas interface and the wall. As a function of polymer fugacity this angle exhibits discontinuities 
of slope ("kinks") upon crossing first-order surface phase transitions located on the gas branch of 
the bulk binodal. Each kink corresponds to a transition from n — 1 to n colloid layers adsorbed at 
the wall, referred to as the n'th layering transition. The corresponding adsorption spinodal points 
from n — 1 to n layers upon reducing the polymer fugacity along the bulk binodal were found in a 
previous study (J. M. Brader et al. J. Phys.: Cond. Matt., 14: LI, 2002; Mol. Phys., 101: 3349, 2003). 
Remarkably, we find desorption spinodal points from nton-1 layers to be absent upon increasing 
polymer fugacity at bulk coexistence, and many branches (containing up to 7 colloid layers) remain 
metastable. Results for the first layering binodal and both spinodal branches off-bulk coexistence 
hint at a topology of the surface phase diagram consistent with these findings. Both the order of 
the transition to complete wetting and whether it is preceded by a finite or an infinite number of 
layering transitions remain open questions. We compare the locations of the first layering binodal 
line and of the second layering binodal point at bulk coexistence with recent computer simulation 
results by Dijkstra and van Roij (Phys. Rev. Lett., 89: 208303, 2002) and discuss our results for the 
contact angle in the light of recent experiments. 



I. INTRODUCTION 



In suspensions of sterically-stabilized colloidal particles 
mixed with nonadsorbing globular polymers, the latter 
induce an effective attraction between the colloids due 
to the depletion effect Q . Such mixtures can phase sep- 
arate into two fluid phases, one being a colloidal liquid 
that is rich in colloids and poor in polymers and the other 
being a colloidal gas that is poor in colloids and rich in 
polymers. Colloid-polymer mixtures serve as excellent 
model systems to study many phenomena associated with 
liquid-gas phase separation as time and length scales are 
much larger than in atomic and molecular systems Q, • 
Recent experiments have focused on the bulk phase be- 
haviour p| (and Refs. therein), (colloidal) liquid-gas in- 
terface tension 0, 0, |Mfl , capillary wave fluctuations 
observed in real space @], droplet coalescence and fur- 
ther non-equilibrium phenomena P,Q,I3- The behaviour 
at nonadsorbing walls has been studied by measuring the 
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contact angle between the free liquid-gas interface and a 
substrate acting as a hard wall. Complete wetting of the 
wall by the colloidal liquid has been observed for a wide 
variety of statepoints However, there are also reports 
of a transition from partial to complete wetting jlOl lllj , 
hence this remains an interesting topic. 

Most of the essential physics of colloid-polymer mix- 
tures is captured by the Asakura-Oosawa-Vrij (AO) 
model of hard-sphere colloids and ideal polymers [12, 
IT^l . , and which has become a widely-used reference 
system. Theoretical approaches ITa . Il6l and computer 
simulations [H E El |U EJ El El] have given in- 
sight into its bulk phase behaviour, and some recent 
work [2£j, 13 ahns at including more realistic colloid- 
polymer and polymer-polymer interactions. Studies 
based on a one-component description of colloids inter- 
acting with an effective depletion potential (obtained by 
integrating out the polymer degrees of freedom and trun- 
cating at the pair-wise level) were devoted to inhomo- 
gen eous situations, such as the free fluid-fluid interface 
|25l |26| and adsorption at a hard wall |27|]. Following 
the development of an accurate density functional the- 
ory (DFT) specific for the binary AO model [H l^ . 
further research was stimulated in inhomogeneous sit- 
uations such as liquid-gas j^S H^] and wall-fluid inter- 
faces [3(], IS an d results were compared to those from 
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simulations |2J, |2J, |23j, |33J . In particular, in Refs. |30|,|31| 
the AO model was considered in contact with a planar 
hard wall. A sequence of first-order layering transitions 
was found on the gas branch of the (liquid-gas) binodal 
upon reducing the polymer fugacity. Further reducing 
the polymer fugacity leads to a transition to complete 
wetting of the wall by colloidal liquid. This scenario was 
corroborated by a simulation study |2l| . The relation of 
the results from these different approaches will be reex- 
amined in the light of the findings of the present study 
in more detail below. The adsorption properties at a 
wall are intimately related to the wall-fluid interfacial 
free energies (or "wall tensions" ) , for which an analytical 
expression was obtained from a scaled-particle treatment 
and which was found to compare well with results from 
full numerical DFT calculations J32. However, despite 
its experimental accessibility 0, ISL llCt Hl| , the contact 
angle of the free liquid-gas interface and a hard planar 
wall has not been considered neither by theory nor sim- 
ulations, in contrast to wetting behavior |2lL l30l l3ll l34| . 
The aim of the present study is to obtain a quantitative 
understanding of the contact angle and elucidate its re- 
lation to the surface phase behavior on the basis of the 
AO model. 

We obtain the (macroscopic) contact angle, 9, from 
the liquid-gas (Ig), the wall-gas (wg) and the wall-liquid 
(wl) interface tensions, 7; s , j wg and j w i, respectively, via 
Young's equation [3{|, 



cos 9 



In 



Iwl 



Hg 
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A similar depletion attraction that acts between two col- 
loids acts between one colloid and a hard wall [27]]; the 
latter therefore favours the colloidal liquid. As a conse- 
quence, we expect j wg > j w i everywhere at coexistence, 
and 9 < ir/2. For statepoints where j wg < j w i + ji g , the 
contact angle 9 > and the surface is partially wet by 
the liquid. However, as soon as j wg = j w i + "fi g a macro- 
scopic liquid layer will intrude between the gas and the 
wall and the latter is completely wet by the liquid. The 
transition from partial to complete wetting induced by 
changing an appropriate thermodynamic variable is re- 
ferred to as the wetting transition |3(| [13 ■ A study of 9 
can supply a link between theoretical predictions of sur- 
face phase behavior and experiments and we display our 
central result in Fig. 

The paper is organized as follows, in Sec. II we define 
the model and and discuss in Sec. Ill the density func- 
tional theory used to calculate the interface tensions. In 
Sec. IV we present results and we conclude in Sec.V. 



II. MODEL 

We consider the Asakura-Oosawa-Vrij (AO) model of 
N c hard-sphere colloids and N p ideal polymers in a vol- 
ume V. The colloids (species c) and polymers (speciesp) 
have diameters <r,, bulk packing fractions r)i = NiVi/V, 



and particle volumes V% = (tt/6)(t| for i = c,p, re- 
spectively. The colloid-colloid as well as the colloid- 
polymer interaction potentials are those of hard spheres, 
so Uij(r) — oo when r < ((7j + <Tj)/2 and Uij(r) = 
otherwise, with ij — cc, cp. The polymers do not in- 
teract with each other, i.e. u pp (r) = for all r. Due 
to the non-additive ranges of these interaction poten- 
tials, the polymers induce an effective attraction between 
the colloids, which for sufficiently large size ratios a v ja c 
(> 0.35) drives a thermodynamically stable phase sepa- 
ration into a colloid-rich (liquid) and a colloid-poor (gas) 
phase [Til fl(| . All bare interactions are of an entropic 
nature, and therefore the temperature T does not play a 
role. The only relevant model parameter is the size ratio 
q = <7p/cr c . We often use the so-called polymer-reservoir 
representation where the mixture is in contact with a 
polymer reservoir, which determines the polymer chemi- 
cal potential. In this situation, the thermodynamic state 
parameters are r\ c and rjp, the latter being the polymer 
packing fraction in the reservoir which is proportional to 
the polymer fugacity as these particles are ideal. 



III. DENSITY FUNCTIONAL THEORY 

We use the fundamental measure density functional 
for the AO model [H, [2^ to calculate colloid and poly- 
mer density profiles from which the interface tensions can 
then be obtained. In density functional theory (DFT), 
the grand-canonical free energy is expressed as a func- 
tional, Q[p c (r), p p (r)], of the one-particle distribution 
functions pi(r) (with i — c,p), given by [3^| 

fi[p c (r),p p (r)} = F cxc [p c (r),p p (r)] 

+ k B T J dv Pl {v) [In (ft(r)AO - 1] 

i=c,p 

+ X] J dr P^ r ) [ u cxt,i(r) - fa] , (2) 

i— c,p 

where fcs is Boltzmann's constant, A, is the "thermal vol- 
ume" of species i, i.e. the third power of the de Brogue 
wave length, and M e xt,i(r) and fa are the external po- 
tential and the chemical potential for species i, respec- 
tively. The excess Helmholtz free energy functional, 
^exc = J dr<f>(r), with $ the excess free energy density, 
is given in Refs. [2£|[2!| and is not reproduced here [47| . 

In thermodynamic equilibrium, the functional is sta- 
tionary, 8Cl/8pi(r) — (with i — c,p), and the resulting 
equations yield the stable distributions, i.e. 



pi(r) = z,exp 



-/3"ext,i(r) - P 



6F cxc [{ Pj (r)}} 



(3) 



with Zi = A" 1 exp[/3/ij] the fugacity of component i and 
f3 = 1/k^T. The equilibrium distribution functions are 
normalized, J drpi(r) — JVf. We note that the polymer 
fugacity is proportional to the polymer packing fraction 
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in the polymer reservoir, rf v = z p V p , and we usually re- 
fer to rf v as the polymer fugacity. We have solved these 
equations numerically for pi{z) in one spatial dimension 
z for the free liquid-gas interface (where u cxt ,i(z) = 
everywhere) as well as both for the liquid and the gas 
at bulk coexistence in the presence of the external hard- 
wall potential, i.e. for u ex t,i(z) = oo for z < (Tj/2 and 
u cx t,i(z) = otherwise for both species i = c,p, where z 
is the space coordinate perpendicular to the wall. The 
numerical routine we have used is a Picard iteration pro- 
cedure with a Broyles mixing scheme |39j. The "mix- 
ing parameter" is continuously adapted to obtain opti- 
mal convergence. Additionally, it is important to realize 
that in situations with several metastable minima, as we 
find to occur for the coexisting gas in contact with the 
hard wall, the initial guess for the profiles in the itera- 
tion procedure determines to which minimum the routine 
converges. 

Once the density profiles are known, the interface ten- 
sion is given by 7 = (^; n h + PV)/A, where Sl; n h = 
n[p c (r), Pp(r)] (i.e. the functional, Eq. [21 evaluated at 
the solutions of Eq. |3J) is the grand-canonical free energy 
of the inhomogeneous system, P is the bulk pressure, and 
A is the lateral (perpendicular to the z-direction) system 
area. In terms of density profiles this quantity can be 
written as 

7 = / dz[w(z)+P], (4) 



where 

<j(z) = k B T ^ Pi( z ) [\n(Pi(z)Ai) - 1] 

i— c,p 

- ViMz) + h B T${z) (5) 

i— c,p 

can be viewed as a local grand potential density (eval- 
uated with the minimized density profiles). In case of 
the liquid-gas interface, the integral in Eq. Q]is over all 
space, i.e. z runs from —00 (bulk gas) to 00 (bulk liq- 
uid). In case of the fluid in contact with a hard wall, 
the integral runs from z = (at the actual location of 
the hard wall) to 00 (bulk). In the numerical routine, 
we compute the interface tensions for each statepoint for 
different system sizes and refined termination criteria for 
the iteration. This gives an estimate of the numerical 
error in the result for the interface tensions, which is 
important as the resulting contact angle can be very sen- 
sitive to these errors. Often the dividing surface |35j is 
chosen at z = <r c /2; then subtraction of J2i=c p ^ ><J i from 
the present definition of the wall-fluid tension, Eq. 0] is 
required. For our present goal this is irrelevant as this 
term does not affect 8, as it drops out of the numerator 
in Eq. □ 



IV. RESULTS 

A. Layered states at the hard wall 

We have calculated colloid and polymer density pro- 
files of the free liquid-gas interface, as well as the wall- 
gas and the wall-liquid interfacial profiles at bulk coex- 
istence. For a given statepoint the liquid-gas and the 
wall-liquid profiles are both unique (disregarding trivial 
translations of the liquid-gas interface); such results have 
been presented elsewhere (3^, El- However, in case of 
the coexisting gas in contact with the hard wall, we have 
found many metastable states, each corresponding to an 
integer number of layers, n, of colloids adsorbed at the 
wall. In Figs.n an d|3 we have plotted a number of such 
profiles, denoted by p c ,n(z)i for size ratio q — 0.6 and 
fugacities r\ r v = 1.2 and %f v = 0.8, respectively, along with 
the corresponding polymer profiles, p v , n (z\ given in the 
respective insets (for reference, the bulk phase diagram 
for q = 0.6 is given in the inset of Fig. |5J. For 77^ = 1.2 
the 0-layer state (marked with an asterisk in Fig. ex- 
hibits practically no excess colloid adsorption and is the 
globally stable state. The grand potential and hence 
the interface tension increases with the number of lay- 
ers, i.e. j wgfi < j W3l i < j wg>2 < ■ ■ • < J wg ,7, where 
is the wall-gas tension corresponding to n colloid layers 
(as given in Eq. 01 and evaluated with Pi. n (z\i = c,p). 
For rfp = 0.8 the equilibrium profile is given by pi_±(z) 
(marked with an asterisk in Fig. Remarkably, in this 
case, the solution pi^{z) corresponds to a higher tension 
than all others, i.e. 7^1 < 7^2 < • • • < Iwgfi < Jwgfi- 
The two state points considered (77^ = 0.8, 1.2) are at 
polymer fugacities larger than that at which the Fisher- 
Widom line liollill hits the bulk binodal (at 77^ w 0.533 
for q = 0.6 |29|), which implies that correlations de- 
cay asymptotically in an oscillatory fashion in the liq- 
uid phase. Apparently, this oscillatory nature also ap- 
pears in the effective interface potential between wall and 
the liquid-gas interface, yielding many metastable min- 
ima H3, see Refs. H2- 

B. Wall contact angle of the liquid-gas interface 

We have calculated the liquid-gas interface tension, 7; g , 
the wall-liquid tension, j w i, and the wall-gas tensions for 
all 71- layer states identified, j wg . ni for the full range of 
polymer fugacities at bulk liquid-gas coexistence. The 
results have been inserted into Young's equation, Eq. ^ 
yielding a contact angle curve, 9 n , for each n-layer state. 
The results for q = 0.6 are plotted in Fig. [21 and a mag- 
nification of the region close to cos 9 n — 1 is displayed in 
Fig. 01 From the definition of the contact angle, Eq. ^ 
and the fact that ji g and j w i are unique for each state- 
point, the state n with the lowest free energy also 
possesses the lowest value of cos#„. Hence, for high 
77p the equilibrium contact angle is given by $o (corre- 
sponding to the 0-layer state), see Fig. 01 Decreasing 77!", 
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FIG. 1: Colloid density profiles, p c .n(z)V c , at a hard wall as a 
function of the scaled distance, z/a c , from the wall at the gas 
branch of liquid-gas binodal for q = 0.6 and rf p = 1.2. Shown 
are results for n-layer states with n = 0, 1, 2, . . . , 7 (left to 
right). All states are metastable except for the globally stable 
n = 0- layer state (see the small peak at contact with the wall, 
marked with an asterisk) . The inset shows the corresponding 
polymer profiles (also from left to right, with asterisk marking 
n = 0- layer state), p P ,n(z)V p , as a function of z/a c . The 
normalizations of the density profiles are such that in bulk 
they reduce to the packing fractions, lim z ^oo pi,n(z)Vi = r\i. 
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FIG. 2: Same as Fig. [I] but for rf p = 0.8 and n = 0, 1, 2, . . . , 6 
(left to right). All profiles are metastable except for the glob- 
ally stable n = 1-state (marked with asterisks). 



leads to an increase in cos#o, until it crosses the cos#i 
branch, and becomes metastable. The crossing point, 
where 0q = 0\, denotes the 0-1 layering transition and 
is also referred to as the first layering transition. (Con- 
sistent with Refs. j^SI^I, the surface phase transition 
from n — 1 to n adsorbed colloid layers at the wall-gas 
interface is referred to as the n'th layering transition.) 
Upon further decreasing r/p (see Fig.0J, cos#i is in turn 
crossed by cos 02 and the crossing point, where 0\ =02, 
is the 1-2 (second) layering transition. This scheme sug- 
gests that there could well be further layering transitions 
upon reducing ry^, and this impression is strengthened 
by the fact that the (metastable) states cos^n (with n 
from 3 to 7, see Fig. |3J|, all seem to converge to around 



the location of the wetting transition. However, we have 
not been able to resolve the third [33, El an d (possible) 
higher layering transitions and we note that any of these 
should be located in the small region of cos0 between 
0.995 and 1 with rf v between 0.6 and 0.65, see Fig. ^ 
Consequently, we can also not obtain insight into the na- 
ture of the transition to complete wetting, i.e. whether 
this is second-order and occurs via an infinite sequence of 
layering transitions or it is first order and is preceded by 
only a finite number of layering transitions. For a more 
extensive discussion of these two possible scenarios, we 
refer the reader to Ref. [3l|. 

Upon reducing 77^, metastable 0-layer states can be 
tracked into a region where the contact angle takes un- 
physical values, cos#o > 1 (Fig. [3J). The same hap- 
pens for 1-layer states as can be seen in Fig. 0] In- 
serting equilibrium values of the interface tensions ob- 
tained from DFT into Eq. ^ ensures that cos 6 < 1, 
but this does not need to be the case when using in- 
terface tensions of metastable states. Reducing 77^ even 
further, the metastable 0-layer state eventually becomes 
unstable at an adsorption spinodal point, see Fig. [3] 
Beyond this point (for even lower ryp, no 0-layer state 
can be stabilized and the numerical iteration rather con- 
verges to the 1-layer state. Similar adsorption spinodal 
points were found for higher n and we have located those 
with a (moderate) resolution of 0.05 in rf p . ^From the 
present data the adsorption spinodal fugacities for n- 
layer states, are rf p<0 , = 0.75, if p V = if p 2 , = 0.65 

and rjp 3 , = jyl 4 , = 77^ 5 , = 0.6. For states with an even 
thicker colloid film (n = 6 and 7), the profiles no longer 
converged properly at low 77^, and we can have not been 
able to obtain precise values for 77^ 6 , and 77^ 7 , . The evo- 
lution for surface states upon increasing rf v is in striking 
contrast. No spinodal points were found and each 71- 
layer state remains metastable up to 77^ = 1.5, a value 
close to the liquid-gas-crystal triple poin t according to 
free- volume-theory for the AO model [lfj . Moreover, for 
large values of 77^ the numerical routine converges very 
rapidly, which hints at deep (nevertheless metastable) 
free-energy minima for these layered surface states. 

As the layering transitions are thermodynamic surface 
phase transitions, they manifest themselves as discontin- 
uous jumps in the Gibbs adsorption 30, 31], 



dz(pi(z) - Pi(oo)), 



(6) 



for both components i = c, p, and these can be obtained 
from Eq. [S] and 



r, = - 
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where "f w f is the wall tension of the fluid. Moving along 
the gas branch of the liquid-gas bulk binodal ties together 
changes in both chemical potentials: 



(8) 
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where the slope of the bulk binodal fulfills a Clapeyron- 
type equation, 



dp c 
dfip 



App 
Ap c ' 



(9) 



which can be deduced from the Gibbs-Duhem equation 
in a straightforward fashion ^3|- Here, Api — p\— p\ is 
the difference in density of species i = c,p in the liquid 
and in the gas phase (note that Ap p < 0). Hence, we 
obtain 



wg 



dp, p 



= -T n 



r, 



App 



(10) 



where the adsorptions T c and T p refer to those of the gas 
at bulk coexistence. Hence, as crossing a layering transi- 
tion (again at bulk coexistence) is necessarily accompa- 
nied by a jump in the adsorptions consequently via 
Eq.^Jthis leads to a discontinuity of slope of the wall-gas 
interface tension. ^From Young's equation it follows 
that this also leads to a jump in the slope of the contact 
angle, d cos 9 jdr\ r v . This is consistent with our findings 
above of crossing of different branches, Q n -\ and 9 n [48| . 

Note that all quantities in Eq.^Jcan be independently 
obtained from our DFT results, namely the adsorption 
Tj from the integral over the respective density profile, 
Eq. EH the differences Api from the bulk phase diagram, 
and the left hand side of Eg. 1101 from a numerical deriva- 
tive of the results for j wg , as obtained through Eq. 0] 
As a check for internal consistency of our calculations 
we have chosen the statepoint with = 1.1 for q = 1 at 
bulk coexistence, which is very close to the binodal of first 
layering transition (which is at rf v « 1.104, as discussed 
below). We find Eq. ^]to be fulfilled to three significant 
digits; for the 0-layer (1-layer) state either sides evaluate 
to 1.264/cr^ (1.352/tr^). Together these two estimates 
yield a jump in dcosO/d-q^ of 0.247, which is consistent 
with our data for the contact angle for q — 1 (which gives 
0.238) and which is discussed below. 

We will now consider the case of larger polymer-colloid 
size ratios, q = 1, see Fig. [5] in more detail. The 0-layer 
branch is calculated for a number of fugacities from high 
values, rf v — 3, to the spinodal point close to ff p = 1. The 
higher-?! branches are only calculated between 0.95 and 
1.5 (n — 1) and between the wetting transition |3l| and 
the crossing point of the 0-layer and the 1-layer branch 
(n = 2 — 6, which practically fall on top of each other). It 
is found that for this size ratio the 0-layer branch is the 
stable everywhere except in a small regime of rf v between 
0.85 and 1.1 (where cos#„ is close to 1) and where layer- 
ing transitions are located. We have been able to identify 
the locations of the first and second layering transitions, 
higher transitions are eroded by the numerical noise. 

Fig. shows the resulting equilibrium values of cos# 
(see the inset for the bare angle 6) as a function of the 
difference in colloid packing fractions of the coexisting 
liquid and gas phases, r) l c — r)9, for both size ratios consid- 
ered, q = 0.6 and q = 1. This representation enables one 
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FIG. 3: Cosines of the contact angle, cos#„ as denned via 
Young's equation as (^ wg —^wi)Hig, as a function of the poly- 
mer reservoir packing fraction, r[ v , for size ratio q = 0.6. A 
scale of the bare angle (in degrees, °) is given on the right ver- 
tical axis. Shown are branches, 9 n , corresponding to n colloid 
layers adsorbed at the wall-gas interface at liquid-gas coex- 
istence (with n = 0, 1, 2, . . . , 7, increasing in the direction of 
the arrow). The constant cos# n = 1 (8 n = 0°) is marked oo 
and corresponds to two macroscopically separated wall-liquid 
and liquid-gas interfaces. For every value of rfp, the branch 
corresponding to the lowest cos6 n (i.e. the largest contact 
angle 8„) yields the thermodynamically stable state; all other 
branches are metastable. The crossing point between 8q and 
9i (marked 0-1) denotes the first layering transition and the 
spinodal point of the 0-layer branch is marked 0'. Also indi- 
cated is the location of the wetting transition (star) according 
to Brader et al. I30II3I. 
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FIG. 4: Same as Fig. El but magnified close to COS Ufi — 1. For 
the sake of clarity, some of the points of the 3-, 4- and 5-layer 
branches close to the wetting transition are omitted. Indi- 
cated is the location of the second layering transition (marked 
1-2) and the 1-layer spinodal point (1'). 



to make direct contact with experiments, where density 
differences of coexisting phases are rather directly acces- 
sible 0, 0, fiol Hl| . For both size ratios the contact angle 
is of the same order of magnitude, i.e. cosO = 0.8 — 1, 
but for q — 1 the wetting transition lies closer (in this 
representation) to the bulk critical point. Only the first 
layering transition for q = 0.6 occurs at a considerably 
large contact angle of cos 6>o-i « 0.95 or #0-1 ~ 18°. A 
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FIG. 5: Same as Fig. |3] but for q = 1. The inset displays a 
magnification of the area close to the wetting transition. The 
points marked 0-1 and 1-2 are the first layering and second 
layering transitions, respectively, and spinodal point of the 
0-layer branch is denoted with 0'. 




FIG. 6: The equilibrium contact angle, cos# (main figure) 
and 8 (in degrees, displayed in the inset), as a function of the 
difference in colloid packing fractions, rf c —r\^^ of the coexisting 
bulk liquid and gas phases for two size ratios, q = 0.6 and 
q — 1. The locations of the equilibrium layering transitions 
are marked by filled circles. The stars are the transitions from 
partial to complete wetting according to Brader et al. [3f3.l3l| . 

remarkable fact is that for both size ratios, there is a sub- 
stantial region in the partial wetting regime close to the 
wetting transition where the contact angle remains very 
small, i.e. in the range (rf c — rf c ) « 0.28 — 0.35 for q = 0.6 
and (r] l c - jg) « 0.2 - 0.27 for q = 1. 

C. Surface phase behaviour 

Layering transitions of the AO model at a hard wall 
were reported earlier by Brader et al. [3(1 l3ll as ob- 
tained within the same DFT approximation |28ll29| and 
using numerical routines similar to those employed in 
the present work. However, in contrast to the present 
study, Brader et al. determined (implicitly) spinodal 
points along two paths: First, by reducing rf v at bulk 



(liquid-gas) coexistence and following the evolution of 
the 0-layer state, observing jumps to 1-layer and sub- 
sequently higher-layer states, and second, by keeping rf v 
fixed and approaching the gas branch of the binodal by 
increasing rj c , starting from very small values. Results 
along reversed paths for the second case, i.e. decreasing 
fi c at constant rf^, were taken to ascertain that hysteresis 
effects are small, and it was concluded that the spinodal 
points give reasonable indications of the locations of the 
equilibrium layering transitions [sil liij . We will discuss 
in the following the relation of these findings to those of 
the present work. For this purpose, we have chosen a ref- 
erence case, i.e. the first layering transition for q = 1, and 
mapped out its surface phase behaviour completely, i.e. 
including binodal and spinodal lines off-bulk coexistence. 
The result is plotted in Fig. [7| 

The layering transition binodal is obtained by deter- 
mining, for each value of 77I, the value of r\ c at which the 
0-layer and 1-layer states have equal interface tensions 
and are thus in thermodynamic coexistence. The result- 
ing binodal line extends below the first layering transition 
at bulk coexistence, at J/po-i ~ •'••■'■) ^° l° wer values of 
rf v into the bulk gas phase region. As this is a transi- 
tion between two (surface) phases of the same symmetry, 
a van-der-Waals loop in the free energy and a critical 
point are mandatory. In order to find the location of 
the critical point, we fit our data for the colloid adsorp- 
tion T c (i]p) of the coexisting 0-layer and 1-layer states 
with a fourth-order polynomial, rf v = ciq + a^C^c — ai) 2 + 
03 (r c — ai) 3 + 04 (r c — ai ) 4 , where the are free fit pa- 
rameters, see Fig. [7| The value of ao is an estimate of 
the critical value of 77^, and the functional form is cho- 
sen to yield the mean-field critical exponent of 1/2, i.e. 
(r c — ai) ~ (r)p — do) 1 / 2 , as the DFT is a mean-field theory 
in the sense that it does not capture fluctuation effects. 
The accuracy in colloid packing fraction r/ c of points on 
the layering binodal is high; typically smaller than 0.1% 
of r/ c . However, the critical point is an extrapolation and 
is therefore much more sensitive to errors, i.e. these may 
be as large as 3% in rf p and 1% in r? c . We have also lo- 
cated the 0-layer spinodal by taking paths at constant rf v 
and increasing values of r/ c monitoring the stability of the 
solution under the iteration procedure. The values of r/ c 
where the 0-layer state becomes unstable and converges 
to the 1-layer state defines the "adsorption spinodal" , 
located at larger values at rj c compared to the layering 
binodal. (The precise location of spinodal points is sub- 
ject to the numerical resolution of the step size in r] c , i.e. 
typically 1% in rj c .) Similarly, we have investigated the 
stability of the 1-layer states upon reducing r\ c at con- 
stant rip. This defines the "desorption spinodal", where 
the 1-layer solutions converge to the 0-layer states and 
which is located at smaller values of r/ c as compared to 
the layering binodal. Upon decreasing 77^, both spinodals 
and the layering binodal end at the surface critical point. 
Indeed even for values of rfp quite above the layering crit- 
ical point the adsorption and desorption spinodals are 
very close and we can confirm the finding of Brader et 
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al. that hysteresis effects at constant rfp are small. 

Next, we discuss the various aspects of this surface 
phase transition in relation to bulk liquid-gas coexistence. 
We first note that no intrinsic difference is observed be- 
tween the layering phase transition in the stable gas re- 
gion and that in the two-phase region where the gas is 
metastable. In the following, we consider the three differ- 
ent surface phase transition lines, i.e. the adsorption spin- 
odal, layering binodal and desorption spinodal, and their 
relation to the bulk binodal. First, the crossing point of 
the adsorption spinodal and the bulk liquid-gas binodal 
denotes the spinodal point terminating the metastabil- 
ity region of the 0-layer state upon reducing rf v at bulk 
coexistence. This corresponds to the adsorption spin- 
odal point of the contact angle (see Fig. O, as located 
previously by Brader et al. (as well as other adsorption 
spinodal points at bulk coexistence for different size ra- 
tios, q = 0.6,0.7,1). Our numerical value of rf p agrees 
well with that of Ref. [3l| • Second, the crossing point 
between layering binodal and bulk binodal is the layer- 
ing transition at bulk coexistence (accompanied by a kink 
in the contact angle, as outlined above). This statepoint 
can be seen as a triple point between the bulk liquid and 
two different surface states of the bulk gas. The location 
of this triple point is quite different from the adsorption 
spinodal point at bulk coexistence. Although hysteresis 
for a path at constant rf v is small, this is not the case for 
the path along bulk coexistence, due to the fact that the 
gas branch of the bulk binodal and the layering (spin- 
odal and binodal) lines have very similar slopes. Hence 
the location of any crossing point is very sensitive to the 
precise location of the individual lines. Third, in striking 
contrast to the previous two cases, the desorption spin- 
odal does not cross the bulk binodal, but remains in the 
one-phase gas region for increasing values of rj^ . We have 
checked this for one additional path at rjL = 2, starting 
with a 1-layer profile at bulk coexistence and decreas- 
ing rj c and indeed found the desorption spinodal point 
in the one-phase gas region. This behaviour is consis- 
tent with the behaviour of the contact angle, which we 
discussed above for q = 0.6 and q = 1, and where ad- 
sorption spinodal points were found upon decreasing j£, 
but no desorption spinodal points were found upon in- 
creasing r]p. Although we have determined this scenario 
only extensively for the first layering transition for q = 1 , 
we believe it to hold more generally for higher layering 
transitions and other size ratios. 

In order to summarize our present results and those of 
Ref. [3l| for the layering transitions of the AO model at 
a hard wall, we draw the surface phase diagram for q = 1 
in Fig. [SJ The results for the n-layer adsorption spinodal 
points for n = 0,1,2,3 at bulk coexistence are taken 
from Ref. |3l|. At bulk coexistence, the first and the sec- 
ond layering binodal points are located at higher fugac- 
ities compared to the corresponding 0-layer and 1-layer 
adsorption spinodal points, respectively. Higher layer- 
ing transitions, corresponding to 2- and 3-layer spinodal 
points of Brader et al. do not emerge from our present 



data. On the other hand, Brader et al. have not found 
layering lines extending into the bulk gas phase for q = 1 
and which we have located for the first layering transi- 
tion. We have not searched for a similar layering line in 
case of the second layering transition. However, from the 
topology of the first layering transition, which we estab- 
lished above, and taking into account the considerable 
separation of the second layering binodal and the 1-layer 
adsorption spinodal points at bulk coexistence, it seems 
plausible that such a layering binodal line also exists for 
the second layering transition. 



Next, we compare our results to those from simulations 
by Dijkstra and van Roij for q = 1 )2l|. They find dif- 
ferent regimes of complete and partial wetting, and first, 
second and third layering transitions located off-bulk co- 
existence. Our results for the first layering binodal line 
extending into the one-phase gas region are in qualita- 
tive agreement with their findings. Previous comparisons 
with DFT results [3l| left a puzzle because a layering line 
was found for q = 0.6, but not for q = 1. The length of 
this first layering binodal line obtained from DFT is con- 
siderably larger (in range of rfp) than that found in simu- 
lations |21|. With hindsight, this deviation seems consis- 
tent with that in other situations, as too small values for 
the critical polymer fugacities compared to simulations 
are also found in bulk (see e.g. 23]) or in confinement 
in planar capillaries [45t 146^ . We have not obtained re- 
sults for the higher-n layering binodal lines, so we can not 
compare the DFT for n > 1 efficiently with simulations. 
Nevertheless, the overall agreement is good and the DFT 
seems to capture all relevant effects: i.e. the layering tran- 
sitions at bulk coexistence, a layering line extending into 
the one-phase region and a wetting transition [3l| . Given 
the approximate nature of the free energy functional and 
the fact that the surface phase transitions are governed 
by tiny free energy differences, this is quite remarkable. 



Coming back to q = 0.6, we summarize our results 
and that of Brader et al. [3(], El in Fig. Similar to the 
case of q = 1, the first and second layering transitions 
at bulk coexistence are located at higher fugacities than 
the corresponding 0-layer and 1-layer adsorption spin- 
odal points, respectively. The effect, however, is more 
dramatic than for q = 1, cf. Fig. Brader et al. have 
also found a 2-layer adsorption spinodal point, indicat- 
ing the existence of a third layering transition, whose 
binodal we have not been able to locate. We have also 
determined the first layering line extending into the bulk 
gas phase and find it to be remarkably long. The location 
of the critical point is determined with the same fit pro- 
cedure as described above for the first layering transition 
for q = 1 . The large separation along the bulk binodal of 
the second layering transition and the 1-layer adsorption 
spinodal point suggests again the presence of a second 
layering line. 
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FIG. 7: The full surface phase diagram of the first layering 
transition (from the 0-layer to the 1-layer state) for q — 1 as 
a function of (bulk) colloid packing fraction, rj c , and polymer 
reservoir packing fraction, rjp. Shown is the layering binodal 
(full line), the 1-layer desorption spinodal (dashed line, left), 
0-layer adsorption spinodal (dashed line, right), and the gas 
branch of the bulk liquid-gas binodal (dotted line). The point 
where the layering binodal crosses the bulk liquid-gas binodal 
is the first layering transition at bulk coexistence (filled circle) . 
The point where the 0-layer adsorption spinodal line hits the 
bulk binodal is the 0-layer spinodal point at bulk coexistence 
(open circle), according to Brader et al. (3l|. Both spinodals 
and the binodal terminate at the layering critical point (filled 
diamond) located off-bulk coexistence in the one-phase gas 
region. Note that the 1-layer desorption spinodal runs entirely 
in the one-phase gas region and hence does not cross the bulk 
binodal. The inset shows the colloid adsorption, r c a^, of 
the coexisting 0-layer and 1-layer states (along the layering 
binodal) as a function of polymer reservoir packing fraction, 
r[ v \ the dividing surface is chosen at z = <r c /2 (lower bound of 
the integral in Eq. |SJ. The symbols are obtained from DFT, 
the full line is the fit (see text) and the arrow at the horizontal 
axis denotes the estimated value of rf v at the layering critical 
point. 



V. CONCLUSION 

In conclusion, we have investigated the contact angle, 
9, of the (colloidal) liquid-gas interface and a hard wall 
using the AO model colloid-polymer mixture and consid- 
ering two different polymer-to-colloid size ratios, q = 0.6 
and q = 1. Our results for 9 are obtained via Young's 
equation from independent numerical DFT calculations 
of the liquid-gas, the wall-gas, and the wall-liquid interfa- 
cial free energies at bulk coexistence. At the planar wall- 
gas interface at bulk coexistence, we identify a range of 
different metastable states each corresponding to a num- 
ber of adsorbed colloid layers at the wall. We argue that 
the globally stable state corresponds to the lowest wall- 
gas interface tension, and therefore possesses the largest 
value of 9. For small density differences of the coexisting 
liquid and gas, i.e. close to the bulk critical point, the wall 
is completely wet by colloidal liquid j^S E3 and hence, 
9 = 0. Moving along bulk coexistence, away from the 



FIG. 8: Summary of the known features of the surface phase 
diagram of the AO model at a hard wall for size ratio q — 1 as 
a function of colloid packing fraction, rj c , and polymer reser- 
voir packing fraction, rjp. Shown are the gas branch of the 
bulk binodal (full line) , equilibrium first (filled circle, marked 
0-1) and second (filled circle, marked 1-2) layering transitions 
at bulk coexistence, as obtained from the present work. The 
first layering critical point (diamond) is connected via the lay- 
ering binodal line (full curve) to the first layering transition at 
bulk coexistence. The spinodal lines are omitted for clarity. 
The 0-, 1-, 2-, 3-layer adsorption spinodal points at bulk co- 
existence (open circles, marked 0', 1', 2' and 3' respectively) 
and the wetting transition (star, marked W) are taken from 
Refs. H3]). The inset shows part of the data together with 
the bulk critical point (large filled square) on a larger scale. 




0.01 0.02 0.03 0.04 0.05 0.06 
'c 

FIG. 9: Same as Fig.[H]but for q = 0.6. The 0-layer adsorption 
spinodal |30L l3l| connecting statepoint 0' with the surface 
critical point (diamond) is omitted for clarity. 



critical point, the wetting transition |3Ct l31| is crossed, 
and 9 becomes non-zero. However, typical values of 9 
remain very small in a considerable range of statepoints 
in the partial wetting regime. Mediated by a sequence 
of first-order layering transitions at the wall-gas inter- 
face, the contact angle grows upon moving further from 
the critical point and reaches typical values up to ~ 35°. 
These layering transitions of the coexisting colloidal gas 
in contact with the wall appear as a discontinuities in 
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the slope of 9 as a function of a thermodynamic control 
parameter, e.g. the polymer reservoir packing fraction or 
the colloid packing fraction difference between both co- 
existing bulk phases. 

Previous difficulties to measure the contact angle accu- 
rately have been overcome by Aarts and Lekkerkerker 
with the use of confocal scanning laser microscopy 111 
a system with q = 0.56. The authors conclude that 9 = 
for all statepoints considered, consistent with their direct 
observation of a prominent colloid wetting film at the in- 
terface of the bulk gas with the wall. They point out 
that actual values of 9 are very sensitive to the precise 
determination of the location of the wall. Large values 
of contact angles as well as the observation of the transi- 
tion to complete wetting have been reported by Wijting 
et al. |10L lll| , using extrapolation of dynamical measure- 
ments (i.e. moving the wall) to zero velocity. However, 
some reservations have been made with respect to the 
latter results The magnitude of the contact angle 
results from subtle differences between the interface ten- 
sions and we do not expect our present results to resolve 
experimental issues. The contact angles which we have 
calculated for the highly idealized AO model in contact 
with a hard wall can therefore only serve as a reference 
case. Important effects due to more realistic polymer- 
polymer interactions |34j . polydispersity and gravity are 
not captured in our present model (see Ref. |33| for a 
discussion) . 

We have also reconsidered the surface phase behavior 
of the AO model colloid-polymer mixture at a hard wall. 
This system is known to exhibit a sequence of first-order 
layering transitions upon following the gas branch of the 
liquid-gas bulk binodal towards the bulk critical point 
(i.e. reducing rf v ). In addition, layering lines extending 
off-bulk coexistence into the one-phase gas region have 
been located. Such a layering transition is characterized 
by a jump in the colloid adsorption at the wall and can be 
identified as the growth of an additional colloid layer at 
the wall-gas interface 0, |3(], El ■ For one specific case, 
being the first layering transition for size ratio q = 1, we 
have determined the layering binodal, which gives the 
equilibrium location of the transition, to high accuracy. 
In addition, we have located the 0-layer adsorption spin- 
odal line, beyond which (for higher rj c at constant rfp) 
the 0-layer state is unstable and the I-layer desorption 
spinodal line, marking the end of stability of the I-layer 
state (i.e. for lower rj c at constant rfp). The layering bin- 



odal and the adsorption and desorption spinodal lines 
end at a critical point, located in the single-phase gas 
region of the bulk phase diagram. The crossing point of 
the layering binodal and the bulk binodal represents a 
triple point between the bulk liquid and the two layered 
states (0 and 1 layers) of the bulk gas which have dif- 
ferent values of the adsorption of both components. We 
find the location of this triple point to differ substantially 
from the (previously identified Hi], El) crossing point of 
the adsorption spinodal and the bulk binodal. Remark- 
ably, we could not find a crossing point of the desorption 
spinodal and the bulk binodal and its absence gives rise 
to continued (meta)stability of the 1-layer state upon in- 
creasing rjp at coexistence. We believe that this scenario 
holds for higher layering transitions and other size ratios. 
We have presented further results for the second layering 
transition for q = 1 as well as for first and second lay- 
ering transitions for q = 0.6. The order of the wetting 
transition and whether it occurs via an infinite or a finite 
sequence of layering transitions remain open questions. 
Whether the occurrence of layering transitions is specific 
to the AO model (see also Ref. |31| for a more extensive 
discussion) or would be present in more realistic descrip- 
tions of colloid-polymer mixtures, is another interesting 
question. However, any experimental attempt to reveal 
such (layering) phase behavior would require exceptional 
accuracy for determining 9 or resolution on the particle 
level for direct observation. 
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